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ABSTRACT 

Discovered in October 2010 by the LINEAR survey, P/2010 TO20 LINEAR- 
Grauer (P/LG) was initially classified as an inert Jupiter Trojan. Subsequent ob- 
servations obtained in October 2011 revealed P/LG to be a Jupiter-family comet. 
P/LG has one of the largest perihelia (q = 5.1 au) and lowest eccentricities (e = 0.09) 
of the known JFCs. We report on observations of P/LG taken on 29 October 2011 and 
numerical simulations of its orbital evolution. Analysis of our data reveals that P /LG 
has a small nucleus (< 3 km in radius) with broadband colours (B — R = 0.99 ± 0.06 
mag, V — R = 0.47 ±0.06 mag) typical of JFCs. We find a model dependent mass-loss 
rate close to 100 kg s _1 , most likely powered by water-ice sublimation. Our numerical 
simulation indicate that the orbit of P/LG is unstable on very short (10 to 100 yr) 
timescales and suggest this object has recently evolved into its current location from 
a more distant, Centaur- type orbit. The orbit, dynamics and activity of P/LG share 
similarities with the well known case of comet 29P/Schwassmann-Wachmann 1. 

Key words: comets: individual (P/2010 TO20 LINEAR-Grauer) — methods: data 
analysis — minor planets, asteroids — solar system: general — techniques: photomet- 
ric — methods: n-body simulations. 



1 INTRODUCTION 

Jupiter family comets (JFCs) originate in the transneptu- 
nian region of the solar system known as the Kuiper belt. 
Kuiper belt objects (KBOs) preserve key information about 
the epoch of planetesimal formation. The low ambient tem- 
perature (~ 40 K) allows KBOs to retain most if not all of 
the ices present in their formation environment. This makes 
them valuable time capsules with which to study an epoch 
long gone. JFCs represent the small end of the steep size 
distribution of KBOs. They are survivors of a dynamically 
intermediate population, the Centaurs, that were neither 
ejected from the solar system nor collided with one of the 
giant planets (Jupiter to Neptune) on their journey from the 
Kuiper belt into the inner solar system. 

The surfaces of JFCs are heavily processed when com- 
pared to KBOs, mainly by sublimation of surface ice. Indeed, 
the optical broadband colours of the two populations dif- 
fer substantially: JFCs have typically solar (neutral) colours 
(B — R ~ 1.3, a — r ~ 0.6), with a spread of about 0.2 mag 



a broad range of optic al colours, from neutral (B — R~ 1) to 
very red (B-R ~ 2.5. |jewitt fc LuvJl200ll : |jewitt|[2002l 'l. The 
neutral surfaces are usually attributed to a fresh ice coating 
while very red surfaces are thought to be the end product 
of irradiation of initially neutral material (1 Thompson et all 
1 19871 : IM oroz et aT][l9 98). The intermediate Centaurs present 
a peculiar, also intermediate distribution of colours which 
includes neutral objects (blue group) and v ery red objects 
(red group) but very few cases in between (|Peixinho et alj 
2003). This Centaur co lour bimodality may be related to 
their dynamical history dMelita fc Licandroll2012f) or simply 



Lamv fc T oth 2009; Sol ontoi et afll2012h while KBOs span 
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to the fact that they are small ( Peixinho et al.l l2012). 

Some objects straddle the line between Centaurs and 
JFCs. They are active, like JFCs, but have Centaur-like or- 
bits. A famous example is 29P/Schwassmann-Wachmann 1 
(hereafter 29P). The orbit of 29P is nearly circular (e = 0.04) 
with perihelion and semimajor axis beyond Jupiter (q = 5.7 
au and a — 6.0 au), and obeys the dynamical definition 
of a Centaur. 29P is constan tly active and displays spo- 
radic photometric outburst s (|Roemerl Il958l : Ijewitj Il99d : 
iTrigo- Rodriguez et all 120081 '). The activity of 29P is mainly 
driven by carbon monoxid e sublimation (|Senav fc Jewittl 
1 1994 ICrovisier et"aH Il995h . Comet 29P has a radius r = 
23 ± 3 km (Yan Fernandez, private comm.). 

In this paper we present a study of comet P/2010 TO20 
LINEAR-Grauer (hereafter P/LG) which is in many re- 
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Figure 1. Osculating orbit of P/LG on 3 October 2010 (near the 
time it was discovered by LINEAR) and on 29 October 2011, when 
the observations reported here were taken. Due to the proximity 
to Jupiter the osculating orbit of P/LG changed significantly in 
a period of only one year. An arrow marks the direction to the 
last perihelion, the orbits of Earth, Mars, Jupiter and Saturn are 
plotted as dashed lines and the axes are labelled in astronomical 
units. 



Table 1. Orbital properties of LINEAR-Grauer on 2011 Oct 29. 



Property 



Value 



Scmimajor axis, a 5.610 au 

Eccentricity, e 0.087 

Inclination, i 2? 628 

Argument of perihelion, w 252? 9 

Longitude of ascending node, f2 44? 03 

Mean anomaly, M 84? 22 

True anomaly, v 94? 19 

Last perihelion passage 2008 Sep 5 

Perihelion distance, q 5.122 au 

Aphelion distance, Q 6.097 au 



spects similar to 29P. Comet P/LG has a low eccentricity 
(e = 0.09), high perihelion (q = 5.1 au) orbit, very close to 
meeting the Centaur definition (see Table [TJ Fig. [1} ■ When 
P/LG was discovered by the Lincoln Near Earth Asteroid 
Research (LINEAR) survey, on 1 October 2010, it was mis- 
classified as a Trojan and its weak activity went unnoticed. 
We present optical broadband measurements of P/LG and 
numerical simulations of its orbital evolution. The former 
show that P/LG is weakly active , at a level comparable to 
the active Centaurs l|Jewittl [20091 ) while the latter indicate 
that P/LG has in all likelihood recently evolved into its cur- 
rent location from a Centaur orbit. The implication is that 
P/LG may offer a glimpse of a relatively small and fresh 
Centaur object. 



Table 2. Observing Geometry. 



Property 


2010 October 03 


2011 October 29 


Solar phase angle (a) 


2.07° 


0.89° 


Heliocentric distance (R) 


5.295 au 


5.603 au 


Geocentric distance (A) 


4.309 au 


4.613 au 


1" at distance A 


3135 km 


3355 km 



2 OBSERVATIONS 

P/LG was observed on 29 October 2011 at the ESO New 
Technology Telescope (NTT) located at the La Silla Ob- 
servatory, Chile. The night was photometric and the see- 
ing varied between 0. 8 and 1.0". At the NTT we used the 
EFO SC2 instrument l)Buzzoni et all Il984l : ISnoderass et all 
2008) which is installed at the f/11 Nasmyth focus and is 
equipped with a LORAL 2048 x 2048 CCD. We used the 2x 2 
binning mode to bring the effective pixel scale to 0.24"/pixel. 
Our observations were taken through Bessel B,V,R filters 
(ESO #639, #641, #642, respectively). 

The images of P /LG were collected in a relatively regu- 
lar fashion, in sets of three consecutive exposures per filter, 
with exposure times of 240 s for B and 120 s for V and 
R. In total, we collected 12 images in the B band, 12 in 
V , and 18 in R. Throughout the observations of P/LG, the 
telescope was set to track the non-sidereal motion of P/LG 
at an approximate rate of —18" /hour in right ascension and 
—6" /hour in declination. 

Bias calibration frames and dithered twilight flats 
through all three filters were collected on the same night 
as the science data. The reduction of the science images, 
consisting of standard bias subtraction and flat fielding, was 
done using the IRAF ccdproc routines. The R band images 
suffered from fringing which was removed using an IRAF 
package optimised for EFOSC2 that was kindly supplied by 
Colin Snodgrass. The P/LG data we re absolutely calibrated 
using observations of lLandoltl (|l992T ) stars taken throughout 
the night. 

P/LG was observed as part of the Pan-STARRS PS1 
all-sky survey very near the time of its discovery by LIN- 
EAR. Located on Haleakala, Maui, the 1.8-m PS1 telescope 
is equipped with a 1.4 gigapixel camera covering 3.2° x 3.2° 
on the sky. The survey repeatedly covers the 37r steradians of 
sky visible from Haleakala. PS1 uses a photometric system 
that approaches the SDSS filter system with the addition of 
a wide (w) band filter which roughl y corresponds to th e com- 
bined pass band of the gri filters (|Tonrv et al.l 12012). The 
PS1 gigapixel camera has a pixel scale 0.25"/pixel. P/LG 
was imaged in four consecutive 45 s exposures taken through 
the w filter on 3 October 2010. All images were processed au- 
tomatically by the Pan-STARRS Image Processing Pipeline. 



3 COMET MORPHOLOGY 

Figure [5] shows R- and B-band summed stacks of P/LG and 
Fig. [3] shows enlarged sections of the region surrounding the 
nucleus of P/LG, to highlight details. A description of how 
the stacks were assembled can be found in H4.ll and H4.3I 
here we focus solely on inspecting the overall appearance of 
the comet. The projected comet tail is clearly visible and 
extends for more than 1' in the direction opposite the or- 
bital motion. Due to the very small solar phase angle, the 
antisolar direction points almost radially away from the ob- 
server. As a result, the projected, antisolar component of the 
cometary tail is tiny, and more easily viewed in the enlarged, 
median-filtered version shown in Fig. [3] The bulk of the tail 
particles trails the comet as a result of Keplerian shear. The 
Laplacian-filtered image (Fig. [3} displays a darker ridge at 
the centre of the tail which originates in the nucleus and is 
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Figure 2. Summed stacks of 9 R-band (left) and 6 B-band (right) images of comet P/LG. The projected tail is more than 1' long and 
points in the direction opposite the orbital motion (— v). Because of the very small solar phase angle (Tabled, the antisolar direction 
lies nearly perpendicularly to the plane of the sky. For this reason, the near-nucleus tail has only a small component in the projected 
antisolar direction (—0; see Fig. [3} which rapidly bends towards the direction opposite the orbital motion. 




Figure 3. Enlarged versions of the near-nucleus tail structure from the R-band stack shown in Fig. [5] The original image (left) is shown 
together with versions processed using Laplacian (centre) and median (right) filters. A trail of larger particles leaving the nucleus and 
trailing the comet's orbital motion is visible in the Laplacian-filtered image. The median-filtered image highlights the near-nucleus part 
of the tail which points in the projected antisolar direction. 



3 October 2010, PS 1, w band 



1.6' x 0.9' # ; 



aligned with the orbital motion. The ridge is probably due 
to large particles leaving the nucleus. Larger particles attain 
lower ejection velocities as they couple more weakly to the 
sublimating gas and hence tend to remain concentrated in 
the plane of the orbit where solar radiation acts to push 
them radially outwards and Keplerian shear forces them to 
trail the nucleus. 

Figure [4] shows a stack of four consecutive iy-band im- 
ages of comet P/LG taken by the PS1 survey. The PS1 im- 
ages are taken only 2 days after discovery by LINEAR on 1 
October 2010. The comet appears active, with a tail that ex- 
tends about 12" on the sky. Here, as in the case of our more 
recent observations, the solar phase angle is small, just over 
2° , probably causing the tail to extend mainly into the plane 
of the sky along the line of sight. As above, the near-nucleus 
portion of the sky-projected tail seems to be antisolar, but 
further from the nucleus the tail seems to trail the nucleus 
along the orbital direction. The three sets of images show 
that P/LG was active at the time of discovery (October 
2010) and remains active a year later (October 2011). 



4 PHOTOMETRY 

We performed aperture photometry of comet P/LG and 
nearby field stars using the IRAF apphot task with four 
main goals: 1) to estimate the size of the comet nucleus, 2) 
to estimate the mass-loss rate from the comet, 3) to mea- 
sure the colour of the nucleus and coma dust, and 4) to 
search for brightness variations in the nucleus and coma over 
time. With those goals in mind, we extracted photometric 
data from the science images using four synthetic apertures 



Figure 4. Stack of four consecutive PanSTARRS ui-band im- 
ages of P/LG taken on 3 October 2010, very near the time of its 
discovery by the LINEAR survey. The comet is clearly active in 
these images. 



defining four regions of interest (see Table [3} • The central 
region, 72-4, 4 pixels (0.96") in radius, was selected to rep- 
resent the flux from the nucleus, but includes an unknown 
fraction of light due to the coma. The intermediate annu- 
lus, 72-4,8, between 4 and 8 pixels (0.96" and 1.92") from 
the nucleus includes contributions from both the nucleus 
point spread function (PSF) and the coma. The annular 
region 72s, 13, located between 8 and 13 pixels (1.92" and 
3.12") from the nucleus, is dominated by the coma. Finally, 
the outer annulus 72i3,2o, located between 13 and 20 pixels 
from the nucleus (3.12" and 4.8"), includes coma flux and 
significant noise contamination. Details of how these partic- 
ular apertures were chosen can be found in £|4. II £14.21 and 
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Table 3. Photometry Regions 



Region 


Region 


Aperture radius 


Aperture radius 


Projected radius 


Magnitudes 


Dominant Source 


Label 


Shape 


[pixels] 


["] 


[km] 


in Region 




K4 


Circle 


r = 4 


<t> = 0.96 


d = 3, 220 


7714 


Nucleus (+ Coma) 


^4,8 


Annulus 


r = 4 to 8 


<j> = 0.96 to 1.92 


d = 3, 220 to 6,440 


m4,s 


Nucleus + Coma 




Annulus 


r = 8 to 13 


4> = 1.92 to 3.12 


d = 6, 440 to 10,465 


7718,13 


Coma 


72-13,20 


Annulus 


r = 13 to 20 


<t> = 3.12 to 4.80 


<2 = 10,465 to 16,100 


77713,20 


Coma 



Table 4. Photometry 



Region 


B band 


V band 


R band 


measured 


[mag] 


[mag] 


[mag] 


77(4 


21.94 ±0.05 


21.42 ±0.05 


20.95 ±0.03 


114,8 


22.31 ±0.08 


21.66 ±0.08 


21.21 ±0.05 


"18,13 


22.56 ±0.16 


21.80 ±0.13 


21.33 ±0.08 


777 1 3,20 


22.72 ±0.30 


22.04 ±0.27 


21.47 ±0.15 



H4.3I Magnitudes measured in these regions are labeled 777.4, 
»7l4,8, wi8,i3, and mi3,20. Magnitudes within the annuli, la- 
beled 7ni n . ou t, were calculated from the magnitudes within 
the inner and outer apertures, rrii n and m out , through 

-0.4m out 0.4m;, 



777in out 



-2.51og 10 (l0- 



10" 



We list the B, V and R apparent magnitudes in each region 
in Table H 



4.1 Nucleus Size 

To estimate the radius of the nucleus of comet P/LG we 
began by generating a summed stack of 9 72-band frames 
that were shifted and aligned at the comet nucleus position. 
The stack, displayed in Fig. O has an equivalent integration 
time of 1080 s and offers an enhanced signal-to-noise ratio 
(SNR) detection of the comet. We used the i?-band images 
because they are sharper (better seeing) and have higher 
SNR than those in B and V . The particular frames that 
were stacked have seeing between 0.8" and 0.9". An analo- 
gous stack that aligns the field stars instead of the comet was 
used to investigate the best choice of aperture for the pho- 
tometry. We found that a 4-pixel radius aperture, containing 
~ 85% of the flux of a point source, offers the best com- 
promise between maximising the SNR within the aperture 
and minimising the contribution from surrounding sources. 
We labeled this region 72.4 and used the magnitude within 
it, ni4, as best representing the flux from the nucleus. In 
the case of comet P/LG we obtain a SNR_r ~ 75 within 
the central region 72.4 and measure an 7?-band magnitude 
m 4 = 20.95 ± 0.03 mag. 

Assuming a spherical cometary nucleus with cross- 
section 7rr 2 , we calculate d the equivalent radius r e using 
the relation l|Russelllll91rJ ): 

-0.4(m4 — m 0) 



Tvr 2 e p R W' - 400 ' = 2.25 x 10 22 n B 2 A 2 10 



(1) 



where pn is the i?-band geometric albedo and /3 is the linear 
phase coefficient of the nucleus. Table[2]lists the values of the 
solar phase angle, a, heliocentric distance, R, and geocentric 
distance, A, of P/LG on 29 October 2011. The apparent red 
magnitude of the Sun at Earth is 771q = —27.1 mag. The 



values of pa and (3 are unknown for P/LG. The uncertainty 
introduced by /3 is negligible when compared to that in pu 
(uncertain by a factor 2 or more) because P/LG was ob- 
served at very low phase angle (a < 1°). We used ft = 0.02 
mag/ ° (|Millis. Ahearn fc Thompsonl [i"98l iMeech fc Jewittl 



11987V) and p R = 0.1 (|Kolokolova et alj|2004h and found an 



equivalent radius r e < 3 km. This figure is an upper limit 
due to the unknown contribution of dust coma to the flux 
within the central aperture. 

The apparent magnitude of P/LG may be converted 
to a nucleus absolute magnitude, m(l, 1,0) (the theoretical 
magnitude of the nucleus when placed at 1 au from the Sun 
and the Earth, and seen at 0° phase angle), using: 



m(l, 1, 0) = m 4 - 5 log 10 (J?A) - /3a. 



(2) 



Substituting the quantities R, A and a from Table [2] and 
again using f3 = 0.02 mag/° we obtain mn(l, 1, 0) = 13.04± 
0.03 mag. The formal error in m(l, 1, 0) is the same as that 
in the apparent red magnitude, but the true uncertainty is 
higher due to the unknown value of /3. We also note that 
this is a lower limit to the absolute magnitude of the P/LG 
nucleus because of the unknown contribution of near-nucleus 
coma to the magnitude 7774. 



4.2 Mass-loss rate 

An order of magnitude estimate of the mass-loss rate from 
P/LG can be calculated by dividing the total dust mass 
present within a coma-dominated annulus surrounding the 
nucl eus by the time it takes the dust to cross the annu- 
lus l|Jewitt fc Luulll989r i. As in 34~H we used the 9-frame, 
7?-band stack to measure the amount of sunlight-reflecting 
dust within the annular region. The annulus was set be- 
tween 8 and 13 pixels from the nucleus photocentre and 
labeled 728, 13 (Tabled . The annulus dimensions were found 
by experimentation (using the profiles of field stars in the 
star-centred stack) to ensure that the flux within the annu- 
lus is dominated by the coma and sufficient to overcome the 
noise due to the sky background. In our stacked data, only 
3% of the flux of a point source is present in the wings of 
the PSF beyond an aperture of 8 pixels. 

The 7?-band magnitude within the annulus 72s, 13 is 
7778,i3 = 20.95 ± 0.03 mag (see Table 0J. The magnitude 
7778,13 can be converted to a total dust cross-section, Ad, 
using Eq.[T]if we use the latter instead of the nucleus cross- 
section 7rr 2 . In this case, pn and f3 are the 7?-band albedo and 
linear phase coefficient of the dust particles, also unknown 
for P/LG. We assumed that these values are the same as 
for the nucleus (/3 = 0.02 mag/° and pn = 0.1) and found 
a total dust cross-section area Ad ~ 2 x 10 7 m 2 . This fig- 
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Figure 5. Radial variation of the B — R colour (thick horizontal 
lines) and uncertainty (grey boxes). The average colours (thick 
horizontal lines) are shown in each of the four concentric regions 
described in Table E] The coma-dominated annuli (apertures > 
8 pixels) appear redder than the central region despite the large 
uncertainties. 
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Figure 6. Lightcurvcs of comet P/LG sampling the central re- 
gion 7?-4 (A1TI4, measured within a 4 pixel or 0.96" radius aper- 
ture), the intermediate region 7£4,8 (Am4,8, from 0.96 to 1.92"), 
and the coma-dominated region "R-8,13 (Am8,i3, from 1.92 to 
3.12"). For each region, we plot relative magnitude lightcurves in 
bands B (blue squares), V (green diamonds) and R (red circles). 
The lightcurves in dilferent bands were aligned simply by median 
subtraction, and the lightcurves in each region (Ani4, Am4^, 
Arras, 13) were shifted vertically for clarity of presentation. 



ure is uncertain by at least a factor of 2 mainly due to the 
uncertainty in the albedo. 

To convert the dust cross-section, Ad, to a dust mass, 
Md, we assume that the dust consists of equal-sized spheres 
with an equivalen t radius, Vd, rep resentative of their true size 
distribution (e.g.. lLi et al.l[201ll ). The mass is then given by 



M d = W3)p d r d A d 



(3) 



where pd is the bulk density of the dust particles. We take an 
equivalent dust particle radius of = (0.1 /im x 1 cm) 1 ' 2 ~ 
32 /im which is based on a power-law dust size distribu- 
tion proportional to (dust grain radius) ~ 3,5 with mi nimum 
and maximum grai n radii of 0.1 pm and 1 cm |jewittj 
120091 ; iLi et al.ll201lT ). We further assume a dust bulk den- 
sity pd = 1000 kg m~ 3 and find a total dust mass within the 
annulus M d ~ 8.6 x 10 5 kg. 

The annulus crossing time, t C ross, is obtained by divid- 



ing the annulus projected width, w c 



4 x 10° m (see 



Table [3]), by the dust velocity, Vd- The v elocity Vd is highly 
uncertain and depends on grain size dCrifo et all [2004) 
but estimates based on macroscopic fragment ej ection from 
17P/Holmes l|Stevenson. Klevna fe Jewittll2010|) and on the 



coma expansion velocitie s of 17P/Holmes (Montalto et al 
2008] ; iHsieh et all |2010| ) and C/Hale-Bopp (|Biver e~ 



2002J) vary from a few 100 m s _i to 1000 m s . We take an 



intermediate Vd = 500 m s 1 and obtain an annulus crossing 
time 



tc 



I'd 



0.8 x 10 4 



(4) 



Finally, the mass loss rate is given by dMd /dt = Md /t CTOSS ~ 
1.1 x 10 2 kg/s. At this rate, a nucleus of radius r e = 3 km 
(see H4.ll) and bulk density p nuc = 1000 kg m~ 3 would have 
a sublimation lifetime t BU b — 3 x 10 4 yr, significantly shorter 
than the dynamical lifetime of P/LG (see q5.2[) . 



4.3 Nucleus and Coma Colour 

We investigated the possibility that the nucleus and coma 
dust of P/LG have different colours by comparing the fluxes 
through different filters within the regions detailed in Table 
[3] To do so, we generated B- and V-band stacks as described 
in £14. II using the frames with best image quality (6 x B- 
band, seeing 1.0" to 1.1"; 6 x V-band, seeing 0.8" to 0.9") 
to obtain single images with equivalent exposure times 1440 
s in B (see Fig. [2]) and 720 s in V. The resulting stacked 
images have SNRs ~ 32 and SNRv ~ 34 within the central 
aperture. 

Figure [5] plots the B — R colour of P/LG versus 
projected distance from the nucleus, and Table 2] lists 
the B, V and R apparent magnitudes in the regions 
of interest. The comet nucleus region, IZ4, has colours 
B - R = 0.99 ± 0.06 mag and V - R = 0.47 ± 0.06 
mag, typical of Jupiter family comet n uclei l|jewittl |2002| ; 
ISnodgrass. Lowrv fc Fitzsimmonsl 120061) . Further from the 
nucleus, as the flux becomes dominated by the coma, the 
B — R colour becomes slightly redder. Despite the large un- 
certainties the effect is noticeable and suggests that the coma 
dust may have a different colour from the nucleus. 



4.4 Lightcurves 

We searched for photometric variability in the P/LG data 
using the individual frames taken in B, V and R over the 
course of the night. Periodic photometric variations from 
the central region can signal a rotating, elongated nucleus, 
while temporal variations in the outer coma regions may 
indicate variable activity from the comet. To look for pho- 
tometric variations we performed photometry within the re- 
gions listed in Table [3] for each frame. As part of our observ- 
ing strategy we acquired sets of three consecutive images per 
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filter. That allowed to take the median magnitude of each 
set and divide the median uncertainty by \/3. 

Figure [6] shows the resulting B, V and R lightcurves in 
regions IZ4, IZi^ and 7?.8,i3- Since we are only interested in 
the relative variability, the lightcurves for each band were 
median subtracted for alignment purposes. The lightcurves 
within each region were then shifted vertically for clarity of 
presentation. 

The small variability seen in B and R within the central 
region TZ4 is probably not significant: the magnitudes rri4 in 
B and R display a positive correlation with atmospheric see- 
ing, at a level that has a probability of only 4% of occurring 
by chance. The large scale variation observed in ms : i3, in 
all three bands, does not correlate with seeing nor airmass 
and is probably real. Although the formal photometric un- 
certainties are large, the variation is regular. However, our 
time base is short rendering the interpretation difficult. 



5 DYNAMICAL PAST AND FUTURE 

Discovered by the LINEAR survey in October 2010, P/LG 
was initially classified as a Jupiter Trojan. The object re- 
tained its Trojan status within the JPL Horizons system 
until October 2011 , when cometary activity was detected 
l|Grauer et alj|201ll ). at which point it was reclassified as a 
JFC. Currently, P/LG has a low inclination (i = 2? 6), nearly 
circular orbit (e = 0.09) with a semimajor axis a = 5.6 au. 
The difference between the distances to perihelion and aphe- 
lion is a mere 1 au, leading to a surface temperature varia- 
tion of AT — 15 K (blackbody in equilibrium with sunlight). 
Such a temperature difference seems too small to cause peri- 
odic activity close to perihelion as is the case for most JFCs. 
This issue is discussed in more detail in [JS] 

To investigate the past and future dynamical evolution 
of P/LG, we numerically integrate its current orbit back- 
wards and forwards in time. For this purpose, we use the N- 
body integration package Mercury version 6.2 (|Chambersl 
Il999h . We integrate the orbits of 371 objects, P/LG plus 370 
clones with normally distributed orbital elements centered 
on P/LG's current orbit and with a 1-a dispersion equal 
to 10 -4 of each parametre (Fig. [7}. The Sun and 8 major 
planets are included as massive bodies and the 371 P/LG 
clones as massless test particles. We use Mercury's hy- 
brid algorithm mode which combines a second order mixed- 
variable symplectic algorithm with a Bulirsch-Stoer integra- 
tor to handle close encounters. We select an initial timestep 
of 8 days and remove clones that go beyond 200 au from the 
Sun at any point during the dynamical evolution. 

We note that our simulations of the orbital evolution of 
P/LG neglect non-gravitational forces. Collimated sublima- 
tion jets accelerate the nucleus and may alter the evolution 
presented here. An upper limit to this effect can be calcu- 
lated assuming that the non-gravitational acceleration, T, is 
due to a single sublimation jet, always aligned tangentially 
along the orbit of the comet, through which all mass loss 
occurs. In that case, the rate of change of the semimajor 
axis, da/dt, can be written as 

da/dt = 2Va 2 T/ (GM Q ) (5) 
where V is the orbital velocity, G is the gravitational con- 
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Figure 7. Initial orbital paramctrcs of 370 clones of P/LG used 
in our numerical simulations. The clones are drawn from Gaussian 
distributions with centers on the actual orbital elements of P/LG 
(star symbol) and standard deviations equal to one thousandth 
of the central value. 

stant, Mq is the mass of the Sun, and the acceleration due 
to the jet is 

T = (dM d /dt) (v d /m nuc ) . (6) 

Substituting dMd/dt — 10 2 kg, v d = 500 m s _1 , m nuc 
1,1/ 10 14 kg, and taking V as being the circular velocity at 
a = 5.6 au we obtain da/dt = 1.4 x 10~ 5 AU yr -1 . In the 
sublimation lifetime (t su b — 3 x 10 4 yr) of comet P/LG its 
orbit can decay by (da/dt) x i su b = 0.42 au. This value is a 
strong upper limit, as it assumes a jet geometry contrived 
to achieve maximum orbital decay. 

5.1 Orbit Stability 

Our simulations show that P/LG is dynamically unstable 
on very short timescales. The 1-a spread of the main orbital 
parametres (semimajor axis, eccentricity, inclination) of the 
cloud of P/LG clone particles expands by a factor e in only 
12 to 14 yr into the past and about 100 yr into the future. 
These times provide a reasonable estimate of the Lyapunov 
timescale for P/LG. 

Due to the strongly chaotic nature of the orbit of P/LG, 
we can expect at best to obtain a statistical assessment of the 
past and future evolution of this comet. The chaotic evolu- 
tion also implies that long term int egrations into the past or 
future are statistically equival ent l|Levison fc Duncanlll994l ; 
iHorner. Evans fc Bailevl 2004). However, current close en- 
counters between P/LG and Jupiter (see H5.3P generate an 
asymmetry in the dynamical evolution making it interesting 
to investigate the short term past and future cases sepa- 
rately. Those same close encounters lead to further strong 
divergence of the clone orbits reinforcing the need for a more 
statistical analysis beyond a few kyr from the current time. 

5.2 Orbit History 

When integrating into the past we find that 50% of the 
clones are removed from the simulation by t — —182 kyr. 
The longest lived clone reaches t = —22.5 Myr. The frac- 
tion of surviving clones is plotted versus time in Fig. [8] (left 
panel). Each step in the figure corresponds to the removal 
of 1 clone, so, e.g. only 3 clones survive to times earlier than 
t = — 10 Myr. The Figure also shows the median semima- 
jor axis and perihelion distance of the surviving clones. The 
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Figure 8. Past orbital evolution of the 371 clones of P/LG. The fraction of surviving clones (left) decreases steeply as their orbits are 
integrated back in time. Half the clones do not survive the orbital integration beyond t = —182 kyr. The longest surviving clone reaches 
—22.5 Myrs. The surviving clones evolve into Centaur /Scattered Kuiper Belt Object orbits with increasing median semimajor axis (black 
points in center plot) and median perihelion distance (black points in plot on the right). The gray boxes contain the central 50% of the 
distribution in each time bin. 
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Figure 9. Heat map of the evolution of P/LG clones during 
the backwards integration described in J5] Darker patches have 
been occupied for a longer time by clones. The a — e space is 
divided into 4 regions, representative of Jupiter family comets: 
a < a s/(l + e )i Centaurs: a,s/(l + e) < a < ajv and e < 0.8, Hal- 
ley comets: Centaurs with e > 0.8, and Transneptunian objects: 
a > ajv, where ag and a/v are the semimajor axes of Saturn 
and Neptune. Dashed lines mark the perihelia of Jupiter, Sat- 
urn, Uranus and Neptune. The clones spend most of the time as 
Centaurs, although pathways exist that transport clones to the 
transneptunian region. 



former increases with time into the past, albeit with large 
scatter (indicated by the gray boxes which contain the cen- 
tral 50% of the points). The latter increases to about q ~ 17 
au and then remains approximately constant with a tail to- 
wards larger q. The past evolution suggests that P/LG has 
recently evolved into its current orbit from a Centaur-type 
orbit. 

Figure [5] shows a heat map of the evolution of the 
P/LG clones in the semimajor axis vs. eccentricity plane. 
Darker/redder regions are more often visited by clones as 
their orbits evolve backwards in time. The plane is divided 
into four regions, taken roughly to represent four types of 
orbits: JFC-type, defined as having a < ag/(l + e), Centaur- 
type, bound by as/(l + e) < a < ajv and e < 0.8, Halley- 
type, similar to Centaur-type except with e > 0.8, and 
transneptunian-type, with a > a at, where as and ajv are the 
semimajor axes of Saturn and Neptune. These regions do 
not aim at accurately representing the dynamical behaviour 
of clone particles; they are designed simply as a means to 
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Figure 10. Fraction of surviving clones within each of the regions 
(orbit types) described in Fig. [9] plotted against time. Bins are 1 
Myr wide. The leftmost three bins (—17 Myr to —20 Myr) contain 
a single clone (see Fig- [HJ • 



classify the type of orbit occupied by the clones as a func- 
tion of time. The Figure shows that clones evolve mainly 
along constant-perihelion lines, as their semimajor axis and 
eccentricity increase. Figure [10] traces the history of clones 
across the a-e plane as a function of time. The Figure shows 
the fraction of surviving clones within each of the orbit-type 
regions described above in bins of 1 Myr, 20 Myrs into the 
past. As already seen in Fig. [8] the surviving clones tend to 
quickly move away from the JFC region into Centaur-type 
orbits. Beyond just a few Myrs into the past the semimajor 
axis of the surviving clones moves even further away from 
the sun, into the transneptunian region. The last surviving 
clone possesses a semimajor axis and eccentricity typical of 
a scattered transneptunian object (a ~ 80 au, e ~ 0.65). 

In Figure [11] we show the history of the Tisserand 
parametre of P/LG clones with respect to Jupiter and Nep- 
tune in our numerical simulation. The Tisserand parameter 
is a useful dynamical quantity that is conserved in the re- 
stricted three body problem. Here, the three bodies are the 
Sun, a planet and P/LG. In the solar system, the Tisserand 
parameter is not exactly conserved, but when calculated for 
a given pair of small body and planet it can be used to quan- 
tify the dynamical influence the latter has on the former. For 
instance, the JFCs have Tisserand parametres 2 < Tj < 3 
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Figure 11. Evolution of the Tisserand parametre of clones. 
Myr-wide bins show the fractions of clones that have Tisserand 
parametre 2 < T < 3 with respect to Jupiter, Neptune, both 
planets, or none. A value 2 < T < 3 with respect to a given 
planet indicates that that planet dominates the orbit of the ob- 
ject. The Figure shows that, as time evolves into the past, Jupiter 
hands control to Neptune as main influence to the clone orbits. 

with respect to Jupiter, to which their orbits are strongly 
dynamically coupled. Main belt asteroids have Tj > 3 and 
are relatively stable with respect to the Jovian gas giant. 
The Tisserand parametre of P/LG, calculated with respect 
to planet P, is given by 

T P = a P /a + 2^(l - e 2 )(a/a P ) cos i (7) 

where ap is the semimajor axis of the planet and a, e and i 
are the semimajor axis, eccentricity and inclination of P/LG. 
As shown in Figure 1111 within the first 1 Myr into the 
past the clones of P/LG are controlled partly by Jupiter 
(2 < Tj < 3), partly by Neptune (2 < T N < 3) and partly 
by neither of those two planets. The latter indicates orbits 
between Jupiter and Neptune. Further into the past, Jupiter 
becomes less important dynamically while Neptune becomes 
more dominant. 

In summary, our simulation results conspire to suggest 
that P/LG may be a recent arrival in the inner solar system, 
possibly from a Centaur or transneptunian type orbit. That 
makes this object an interesting target of study as it may 
offer a glimpse of a relatively unprocessed Centaur. 

5.3 Close Encounters with Jupiter and the 
Trojan Clouds 

We performed a high time-resolution simulation (0.1 yr ini- 
tial timestep) going back only 3000 yr to investigate recent 
close encounters between the P/LG clones and the major 
planets. We focus here on Jupiter, due to its current close 
proximity to P/LG. In the very near past, we find that of 
the 371 P/LG clones, 331 (89%) had close encounters with 
Jupiter (< 100 planet radii) between 48 and 51 yr ago. Dur- 
ing that time interval, 34 clone particles came within 20 
Jupiter radii (rj) of the giant planet and the closest ap- 
proach distance for any particle was 13 r.j. Two of the 331 
clone particles had very close approaches to Jupiter in con- 
secutive orbits: one particle came within 14.0 rj just 48 yr 



ago and within 8.7 rj during the previous orbit (60 yr ago); 
another particle passed at 16.0 rj and 14.6 rj about 48 and 
57 yr ago. Finally, if we consider the full 3000 yr time span 
of the simulation, 7 particles came within the Roche limit 
of Jupiter (ss 2.7r\/ assuming a flui d- like P/LG nucleus wit h 
bulk density p^c = 1000 kg m~ 3 ; lAsphaug fc Benz|[l996r i. 
In a similar integration to the one just described but into 
the future we find that 11 out of the 371 clone particles come 
within 10 rj and 3 of those come within 2 rj roughly 145 
yr into the future. A total of 41 clone particles approach 
Jupiter within 100 r.j at that time. 

We also searched for possible crossings of one of the Jo- 
vian Trojan clouds by P/LG clones. The Trojan clouds have 
an irregular shape, roughly 4 x 2 x 0.5 au (FWHM) in siz e 
l|Jewitt. Truiilro fc Luufl2000l ; iNakamura fc Yoshidal 120081 ) . 
To identify potential crossings and estimate the time spent 
by clones inside the Trojan clouds we used Jupiter's posi- 
tion to calculate the positions of its L3 and L4 Lagrangian 
points in each timestep. Finally, for each clone we added the 
timesteps it spent within 0.3 au of either L3 or L4. We chose 
0.3 au to ensure the clone is well within the densest region 
of the Trojan clouds. We found that clones spend an average 
of 10.1 of the last 3000 yr within the regions just described. 



6 DISCUSSION 

Figure |4] shows P/LG as seen by the PanSTARRS PS1 sur- 
vey (see fj2} close to the time the comet was discovered by 
LINEAR. The PS1 images show that P/LG was active at 
the time of discovery and yet the LINEAR pipeline first 
classified this object as a Jovian Trojan. This suggests that 
LINEAR may have misclassified several active objects as in- 
ert asteroids, which is an important point to consider at a 
time when t he study of a ctive asteroids is receiving increas- 
ing interest l|jewittll2012l . and references therein). 

Our data reveal that P/LG was active in October 2010, 
roughly 60° past perihelion, and remains active a year later, 
at 90° past perihelion. P/LG's nearly circular orbit, atypi- 
cal for a JFC, leads us to consider a number of explanations 
for its activity. One possibility is that P/LG has been acti- 
vated by a recent collision with a smaller object. We find this 
unlikely: our simulations show that P/LG spent only a neg- 
ligible fraction of its recent past within the Trojan clouds 
(but see below), and even there the chance of a collision 
onto its r e < 3 km nucleus would be low (intrinsic colli- 
sion p robability P; = 6.5 x 10" 18 km" 2 yr" 1 : Idell'Oro et all 
1998). Another possibility is that P/LG was exposed to sig- 
nificant tidal stress due to close approaches to Jupiter. This 
scenario would require encounters with Jupiter closer than 
the Roche limit (mass shedding begins at d < 0.69J? roc h e 
for a comet on a parabolic orbit; ISridhar fc Tremaindll992l : 
Asphaug fc Bend 1l996l ) which our simulations show to be 
improbable. It is also possible that P/LG became active due 
to rotational mass shedding but since the rotation period of 
P/LG is not known this possibility is untestable. One sce- 
nario that our simulations can not rule out is that P/LG 
originated in the Trojan population and was disloged from 
the 1:1 resonance with Jupiter through the action of non- 
gravitational effects (see SfoJ. Conceivably, P/LG could have 
been collisionally actived by a smaller Trojan and led to drift 
from the stable region, accelerated by mass-loss jets; we note 
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Table 5. Comparison between P/LG and 29P 



Comet 


a 


e 


i 


B-R 


V -R 




M 


dM/dt 


dm/dt 




[au] 




[°] 


[mag] 


[mag] 


[km] 


[kg] 


[kg s- 1 ] 


[kg m — 2 s" 1 ] 


P/LG 


5.6 


0.09 


2.6 


0.99 ± 0.06 


0.47 ±0.06 


< 3 


1.1 x 10 14 


1.1 x 10 2 


7.7 X 10~ 7 


29P 


6.0 


0.04 


9.4 


1.28 ±0.04 


0.50 ±0.03 


23 ±3 


5.1 x 10 16 


5.1 x 10 3 


9.7 x 10~ 7 



Columns are (1) comet ID, (2) semimajor axis, (3) eccentricity and (4) inclination, (5) and (6) broadband 
colours, (7) equivalent radius, (8) model dependent mass-loss rate, (9) specific mass-loss rate. Mass assumes 
nucleus/dust bulk density p n uc — 10 00 kg m" 3 . Specific mass-loss rate assumes a spherical nucleus of radius 
r e . Comet 29P colours from Ijewitj l|2009h . 



that P/LG's highly unstable orbit implies a recent departure 
from the Trojan region. Modelling this possibility is a com- 
plex task involving a number of unknown parametres and 
is beyond the scope of this paper but future observations of 
this comet may shed light on a possible link with the Trojan 
population. 

The activity of P/LG could simply be due to ice subli- 
mation. At first sight, the small difference between the P/LG 
perihelion and aphelion distances (Table [TJ, and the conse- 
quently small orbital variation of the subsolar equilibrium 
temperature (AT « 15 K) could seem insufficient to power 
periodic sublimation activity. However, between the perihe- 
lion and aphelion distances of P/LG, the specific mass loss 
rate due to water ice sublimation varies by more than an 
order of magnitude (from 10 -5 to 10 -4 kg m -2 s _1 ) while 
CO/CO2 sublimation driven mass lo ss varies by a factor of 
2 and is around 10" 2 kg m -2 s" 1 (jjewittj [2009h . We find 
that P/LG loses mass into the coma at a rate ~ 10 2 kg s _1 
implying that an area of 10 6 to 10 7 m 2 would need to be 
active on the surface of P/LG if water ice sublimation is the 
source of the activity. Those areas are small compared to 
the maximum surface area of the nucleus, A-nr\ ~ 10 8 m 2 . 
Sublimation due to CO/CO2 would require only an area 
~ 10 4 m 2 to be active, corresponding to a tiny active vent 
on the surface of the 3 km nucleus. We conclude that given 
the recent perihelion passage and the particular range of 
heliocentric distances traversed by P/LG since then our es- 
timated current mass loss rate is consistent with the activity 
being driven by ice sublimation. 

P/LG is in many ways comparable to comet 29P (see 
Table [S] and Fig. I12p . Both comets have high perihelion 
(q > 5 au), low inclination (i < 10°), nearly circular or- 
bits (e < 0.1) more typical of Centaurs than JFCs. Both 
orbits are unstable on very short tim escales (a few hundred 
years; iHorner. Evans fc Bailevll2004 ). The two comets show 
similar levels of activity, with specific mass-loss rates close to 
dm/dt = 10 -6 kg m -2 s~ , despite 29P being nearly an or- 
der of magnitude larger (>2 orders of magnitude more mas- 
sive) than P/LG. Both objects display continued activity, 
although the coverage of P/LG reported here is not enough 
to establish a pattern of activity. It will be interesting to 
see if P/LG's activity is significantly modulated by its or- 
bital motion or if it remains constant as is the case for comet 
29P. It will also of interest to monitor P/LG and see if it too 
will display the sporadic outbursts seen in 29P. Significant 
differences include the size of the comets and the source of 
activity: 29P is dominated by CO sublimation, while P/LG 
is more likely active due to water-ice sublimation. 
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Figure 12. Comparison between comets P/LG and 29P. Mass- 
loss rates are plotted against nucleus radius (and mass, assuming 
nucleus/dust bulk density p nuc = 1000 kg m — 3 ). The arrow indi- 
cates that we possess only an upper limit on the nucleus radius of 
P/LG. The uncertainty in the radius of 29P is marked by a hori- 
zontal error bar. The size of the P/LG and 29P points is linearly 
proportional to their geometric cross-section for visual compari- 
son. Lines of constant specific mass-loss rate are labeled in units 
of kg m — 2 s — 1 . 



7 CONCLUSIONS 

We report photometric observations and numerical simula- 
tions of the orbital evolution of the unusual comet P/2010 
TO20 LINEAR-Grauer (P/LG). Our findings can be sum- 
marised as follows: 

(i) Comet P/LG was active at the time of discovery (Oc- 
tober 2010) by LINEAR, and remains active in October 
2011. LINEAR did not detect the activity and initially mis- 
classified P/LG as a Trojan, which suggests that several ac- 
tive objects may have gone undetected by the survey. 

(ii) The nucleus of P/LG has equivalent radius r e < 3 
km, and colours B — R — 0.99 ± 0.06 mag and V — R — 
0.47 ± 0.06 mag, values typical of Jupiter family comets. 
The data suggest a slight reddening of the dust colour with 
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distance from the nucleus but the uncertainties are large. We 
find no significant rotational photometric variability from 
the nucleus region. 

(iii) We obtain a model-dependent estimate of the mass- 
loss rate from P/LG of ~100 kg s _1 . We favour water-ice 
sublimation as the simplest and most likely cause for activity 
in comet P/LG. 

(iv) Our numerical simulations show that the orbit of 
P/LG is unstable on very short timescales and suggest that 
it may be a Centaur that recently arrived in the inner so- 
lar system, although other possibilities exist involving non- 
gravitational effects. 

(v) Comet P/LG is in a number of ways reminiscent of 
the well-known 29P/Schwassmann-Wachmann 1. 29P is a 
comet /active Centaur that shows sporadic outbursts super- 
imposed on a background of constant activity. Comet P /LG 
is an order of magnitude smaller (three orders of magnitude 
less massive) than 29P and yet displays similar activity per 
unit area. Comets 29P and P/LG are interesting as possible 
examples of relatively unprocessed Centaurs. 
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